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We show how computational symbolic packages such as FeynArts, FormCalc, Form and 
LoopTools can be adopted for the evaluation of one-loop hadronic electroweak radiative 
corrections for electron-nucleon scattering and applied to calculations of the nucleon weak 
charge. Several numerical results are listed, and found to be in good agreement with the 
current experimental data. 

I. INTRODUCTION 

Next-to-Leading-Order (NLO) effects in electroweak interactions play a crucial role in tests of 
the Standard Model, and require careful theoretical evaluation. An excellent place to search for new 
physics, the deviation of the weak charge of the nucleon from its Standard Model prediction also 
requires considerable experimental and theoretical input. The importance of theoretical predictions 
for the weak charge of the nucleon was recognized more that two decades ago by Marciano and 
Sirlin (3). Their original analysis was followed recently by significant theoretical work on the proton 
weak charge which presented an updated Standard Model prediction for the weak charge of the 
proton, and estimated the QCD structure uncertainties. 

Computer packages such as FeynArts [3], FormCalc 0|, LoopTools 0|, [f| and Form [3] have 
created the possibility to go one step further. Now, we have an option of calculating parity- violating 
NLO effects including all of the possible loop contributions within a given model. We have already 
tested the automated approach to calculate one-quark radiative corrections 0]. 

In the work presented here, we adopt FeynArts and FormCalc for the NLO symbolic calculations 
of amplitude or differential cross section in electron-nucleon scattering. Using Dirac and Pauli-type 
couplings with the monopole form factor approximation, we construct the computational model 
enabling FeynArts and FormCalc to deal with electron-nucleon scattering up to the NLO level. 
Since the weak charge of the nucleon is directly related to the form factors of the parity- violating 
part of the amplitude at the zero momentum transfer, we choose the calculation of the weak charge 
of the nucleon to be a suitable test our method. Within the uncertainty, which mostly comes from 
the uncertainty of the current electromagnetic form factor measurements, our results agree with 
experiment. 

We start with expressions for Pauli and Dirac couplings in terms of fermion weak and electric 
charges, and the definition and classification of one-loop radiative corrections. After that, an 
example for the {7 — Z} box type of correction for electron-proton scattering is considered and a 
computational model is proposed. We proceed with computational details for the self-energy graphs 
and vertex correction graphs. Since we reserve full kinematic dependence in all types of our NLO 
calculations, it will make it easier to adopt our results to the current and future parity-violating 
experiments. 
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II. THEORY 



A. Dirac and Pauli Coupling 



In the approximation where the nucleon behaves as a point-like particle, vector boson couplings 
obey general rules of electroweak theory. For left and right handed fermions, we can use the following 
structure for the {Z — N} type couplings: 



^Z-N 



ie 



Z—N Z—N 



(1) 



where w± = are chirality projectors and represents the coupling strength for the left 

and right handed fermions, respectively. Substitution of w± into Eq.(pQ) yields the vector and axial 
vector representation of the coupling T^_ N : 



ie 



9t N 



7^75 



(2) 



Couplings of vector bosons to fermions derived from the neutral current part of the electroweak 
Lagrangian predict tree-level values of g L and g R given by: 



9l 



Z-N . Z-N i n -2n 

+ g R 1 C 3 - 4s£,Q 



9 Z R - N -9t N 



C ■),; So. 



(3) 
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Here, C3 and Q are twice the value of the isospin and electric charge respectively, and s w and c w 
are sin^ arid cos#vKj where Oyyis a Weinberg mixing angle. In the case when a photon couples to 



the nucleon, g\ N 



7-JV 
9r 



Q and N = ieQj^. 



To accommodate nucleon structure, we have the couplings preserve their vector and vector-axial 
structure, but with the charges replaced by the corresponding form factors. The electromagnetic 
TJ l ~ N coupling has two vector components responsible for static electric and magnetic interactions: 



IT " (9) = ie 



(4) 



where Fx (q) and F2 (q) are the Dirac and Pauli form factors, respectively, and q a is the four- 
momentum transferred to the nucleon. As for T^~ N (q), we have: 



^Z-N 



(q) = ie 



h (q) 7^ + 2^T a,xaqa ^ 2 te) + 91 ( ? ) > 75 



(5) 



with fx (q), fi (q) and gx (q) as weak Dirac, Pauli, and axial- vector form factors. According to the 
first line of the Eq.j3]), form factors fx (q) and fi (q) are expressed as: 



h,2 (q) 



4:C W S 



(6) 



with Fi 



V( P ) 



1,2 



? V(n) 
1,2 



Ff 2 — Fx 2- For gx (q) we have 



9i (q) 



-9 a (q) 



(7) 
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where g p A (q) = —g\ (q) = gA (q) is an axial form factor. To considerably simplify analytical 
expressions, in our computational model we use the monopole structure for form factors expressed 
as 

/P Xf s A 2 {F h2 ,g A }(0) 

{Fi, 2 , 9a} (q) = A2 ' _ g2 , (8) 

which is a reasonable approximation in our case. The value of the parameter A 2 = 0.83mj^ is found 
after the fit of the electromagnetic form factors by monopole approximation in the low momentum 
transfer region. Further details on how couplings where implemented in the computational model 
are given in the appendix. 

B. Definition of NLO Contribution 

In order to define the Next-to-Leading-Order hadronic corrections, we are using the electron- 
nucleon parity violating Hamiltonian in the form proposed by [H]: 

H PV = —j= [Ci N {Uel^lbUe) ^NI^n) + C 2N (u e l^U e ) {unI^Un)] ■ (9) 

Form factors C\n and C 2 n represent perturbative expansion resulting in 

C{l,2}iV = X! C \l,2}N = C {1,2}N + C{1,2}JV + (° 3 ) • ( 10 ) 

i 

The superscript in CK 2 y N represents the order of the perturbation ("zero"- tree level, "one" - one 
loop level (NLO) and so on). Here Ch 2 y N can be defined as an one-loop contribution to the parity- 
violating form factor, normalized by the Fermi constant, Gf = t^t 2 • A calculation of NLO 
form factors is related to the calculations of loop integrals represented by three topological classes: 
box, self energy, and vertex (triangle) graphs. To preserve gauge invariance we have to include all 
the possible bosons of the Standard Model in these topological classes. Taking into account that in 
the t'Hooft-Feynman gauge the contribution coming from the Higgs scalar and gauge fixing fields 
is negligible, we choose to consider boxes and nucleon vertex corrections (triangles) with 7, Z and 
vector bosons only. For the rest of the graphs - self-energies and lepton vertex corrections - 
we accounted for all the possible particles in the Standard Model. Accordingly, we consider NLO 
corrections for every class. 

C. Box Diagrams in e — N Scattering 

The precise formulae for an entire set of graphs are cumbersome, and it is not feasible to show 
them in the present work. Here, we provide some details for {7 — Z} boxes only. Our complete 
calculations are shown in [<J. The full set of diagrams applicable for this case can be found in [Io| . 

According to the Feynman rules (see, for example [ll||), the amplitude for a {7 — Z} box (see 
Fig. (OP) can be written as: 



16tt 4 J \ e " (pa - k! - k 2 + qf - ml 

(11) 

m N + jq+ fa ^_iv \ (9 ua 9 W \ 



{p 2 + q) - m N ) \ Q {p 2 + q - k 2 ) —m 
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Figure 1: 7 — Z box diagram in the electron-nucleon scattering. 



In order to integrate the amplitude in Eq. lfTTI) and the rest of the NLO loop integrals using 
FeynArts and FormCalc, it is necessary first to create a model file with declarations of all the 
fields participating in interaction and second to describe all the possible couplings of the fields on 
both generic and particle levels. For the case of electron-nucleon scattering, the only new particles 
in the model besides the particles of the Standard Model are neutron and proton, so it is rather 
straightforward to declare them in the FeynArts model file. On other hand, the description of the 
couplings according to the Eq.(JIJ and Eq.j5|) in the FeynArts represents a certain problem. Because 
we use the monopole form factor approximation (see Eq.([8|)), it introduces momentum dependence 
in the coupling's denominator. For FeynArts to deal with momentum-dependent couplings with 
momentum dependence introduced in the denominator, these couplings have to be described as 
propagators. Of course, a propagator introduced in the coupling conflicts with the declaration of 
the fields in the FeynArts model files which use the same propagator notation. We can solve the 
problem by transferring the monopole form factor from the coupling part of model to the part 
where fields are described using propagator notation. For boxes that can be achieved by starting 
with the general definition of a four-point tensor integral of rank k in the form 



[{P2 -ki-k 2 + qf - mfj {{jp 2 + q) 2 - mfy) 



(12) 



1 1 

9 2 (P2 + q ~ k 2 f - m\ 



By adding the monopole form factor approximation into the above definition, we obtain a six- 
point tensor integral which can be reduced into a combination of four-point integrals by using the 
following expansion: 



A 4 



D X D 2 q 2 + q - k 2 f - m% (q* - A 2 ) ((p 2 + q ~ k 2 f - A 2 
A 2 



(13) 



(A 2 -m|) V<7 2 -A 2 q 2j y^ p2 + q _ k2 f_ K 2 ( p2 + q _ k2 ) 





Figure 2: Expansion of the 7 — Z box in terms of £7, bZ particles 

Here D ^ D represents the last two terms of the product in Eq. (fT2l) multiplied by the two monopole 
form factors of the box diagram. Explicitly, the right-hand side of Eq. (fl~3j) can be written in the 
form 

1 A 2 1 1 -A 2 1 1 



D X D 2 (A 2 - ml) (p 2 + q - k 2 f - m| (A 2 - m|) <? 2 - A 2 ( P2 + q _ fc 2 ) 2 - m 2 



z 



-A 2 1 1 A 2 

+ 



(tf-m\)q\ P2 + q -k 2 ) 2 -h? (A 2 -m|)g 2 -A 2 ( P2 + 9 _A :2 ) 2 -A 2 



(14) 



Substitution of the last equation into amplitude Eq.JTTJ results in an expansion where each term 
of the sum represents an amplitude constructed for electron-nucleon scattering with the point- 
like nucleon. Also, in this consideration, the couplings between nucleon and vector bosons are 

multiplied by the factor an d the structure of the second, third and fourth terms of 

Eq. (fT4l) suggests the introduction of fictitious vector boson particles {5^,5Z} with fixed masses 
m {8'y,8Z} = A. These bosons have no physical meaning, of course, but they allow us to remove 
momentum dependence in the denominator of the coupling. Now we can model couplings between 
nucleon and vector fields using Eq.(Jl| and Eq.([5|) with a monopole form factor moved into the 
definition of propagators for the new vector fields such as {5~/,5Z} where sz y = q i_ A2 ■ A 
diagrammatic representation of the proposed expansion is given by the set of Feynman graphs in 
Fig.©. 

As for the boxes, we have all we need to complete the automated construction of amplitudes 
in the FeynArts and the calculations in the FormCalc. We should note that Eq. (fTT]) is written for 
the {7 — Z} box diagram only. To have a complete analysis, it is necessary to consider {Z — 7}, 
{7 — Z} and {Z — 7} crossed boxes as well. {Z — Z} and {W — W} box diagrams should be 
considered, too. Our calculations include 36 boxes in total. 



D. Self-Energy Graphs 

In total, 116 self-energy graphs and 6 counterterms contribute to the PV {e — N} amplitude. 
This includes gauge and the gauge-fixing fields, the Higgs field, and virtual leptonic and quark 
pairs in creation-annihilation processes in the loops. Moreover, the vertex {N — V — N} does not 
belong to the loop integrals and plays the role of a multiplicative factor proportional to the coupling 
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Figure 3: Expansion of the self-energy graphs in terms of £7, SZ partners. Here, factor ±J ^ A2 A ^2 ^ 
for the couplings of the box diagrams should be replaced by one. 



used 




Figure 4: Truncated self-energy graph in hadronic vacuum polarization 

defined in Eq.()4]) and Eq.(|5j). In FormCalc the self-energy loop integrals can be evaluated using 
the expansion given in the Fig. ([3]), and then tensor decomposition and tensor reduction techniques 
applied, leaving the final result as a combination of one and two point scalar integrals which we 
compute using Gauss integration subroutines of LoopTools. 

To cancel ultraviolet divergences we employ an on-shell renormalization scheme according to 
We have to assume that quarks in the self-energy loops are free, but it places certain constraints. 
Since quarks are confined, a QCD strong quark-quark interaction should always be considered. It is 
possible to bypass these complications by replacing quarks with pions, or use "free" quarks but with 
adjusted effective masses. Here, we use the second approach with the effective mass of the quarks 
coming from a fit of hadronic vacuum polarization to the measurements of QED cross section of 
the process e + e~ -^hadrons normalized to the QED e + e~ — > fj> + n~ cross section. The real part of 
the renormalized hadronic vacuum polarization satisfies the dispersion relation: 



A a 



n 



hadr ! 



= -Re U had ( s ) = —s ; 

3vr J 4m 2 s' (s' - 



(15) 



with 



hadrons) 



(16) 



being a very well known experimental quantity and used here as an input parameter. 
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Hadronic vacuum polarization Ti-1 ad (s) is related to the truncated {7 — 7} renormalized self- 
energy (See Fig.(j4j) by the following expression: 

^ fer J^ = ^la d (s) + iIrn(^ ferm (s)), (17) 

which can be easily evaluated by employing the free quark approximation. An updated value of 
the dispersion integral, along with a logarithmic parametrization, can be taken from [jjj. A new 
reported value coming from the light quark contribution at s = m 2 z is Aa^ r {fn z ) = —0.02761. 
This value can be reproduced by Eq.JTTJ using the following masses of the light quarks: m u = 
irid = 53MeV (corresponds to A«^ 5 J r {m z ) = —0.027609). Clearly, the values of the light quark 

masses at low-Q scattering processes should be adjusted by using this approach, but with Aa2fJ r (s) 
calculated in the region of ^/s < 4.0 GeV. The simple logarithmic parametrization can be used here 
to extract quark masses at low momentum transfer: 

^< a %(s) = A + Bln(l + C- S ), (18) 

with A, B and C parameters taken from fl3j ]. For low-momentum transfer experiments, the c.m.s. 
energy is \/s < 4.0 GeV, which gives, m u = ma ~ 45MeV. 



E. Vertex Corrections Graphs 

The vertex correction contributions can be split into two classes. In the first class, where the 
electron vertex is at the one-loop level, the amplitude is calculated according to As in the case 
of the self-energy graphs, the hadronic vertex does not belong to the loop integrals, and therefore 
the PV amplitude is constructed according to the SM Feynman rules employed in the FeynArts 
package. Generally, the electron vertex corrections have an infrared divergence at q — > and are 
treated by the bremsstrahlung contribution considered in 0]. 

The second class of the triangle graphs are hadronic vertex corrections. In this case, to evaluate 
the amplitude automatically using the packages FeynArts and FormCalc, we construct an expansion 
similar to that given for the box diagrams. To work out the set of rules for the triangle topology, 
it is sufficient to consider the example in Fig.j5]). 



s 



N N 




s,z 




Figure 6: Hadronic vertex expansion in the terms of Z, 8\Z and 82Z partners. 

For the graph in Fig.©, the amplitude denominator has the structure 
1111 



D X D 2 D Z D 4 (k A -qf-m 2 N (k 2 - q f- m 2 N (k 4 -k 2 y 



rn 



(19) 



A 2 



A 4 



"4 (&4 - k 2 ) - A 2 (q 2 -A 2 Y 



which can be easily expanded into 
1 1 



D\D 2 D^D 4 (k 4 — q) — m? N (k 2 — q) — m 



N 



D 



Z-N 



+ 



B 



8Z-N 



(k 4 -k 2 y-m 2 z (k 4 - k 2 ) - A 2 



(20) 



lim 



C 



Z-N 



(c^ z - N ) 2 (c 5 



{AiM}->A \ q 2 -m 2 z q 2 - Af 
Here, the coefficients B Z ~ N and B SZ ~ N are equal to db ' A ~ 



2 z-n\ 2 ^ 

) 



2 -A 2 



(A 2 — m|) 

and C" 52 ^ 2 can be calculated using the following formulae: 



respectively. C Z ~ N , C 5lZ - Al , 



(c z ~ N ) 



A' 1 



rat 



A?) (m| - A 2 ) 



,Z-N 



A 1 



A? A 2 -A 2 ' 



(21) 



A 4 



rrir 



A? 



A 2 A? -A 2 /' 



The expansion of the amplitude denominator in Eq. pOj) has a simple graphical representation (see 
Fig.©) and suggests, in this particular case of the triangle topology, that a set of virtual particles 
8\Z and 5 2 Z should be introduced in the NLO hadronic vertex corrections. 
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Although ultraviolet divergences are absent in the hadronic vertex corrections due to the addi- 
tional terms proportional to q -J^_^ in the coupling, it is still necessary to compute the wave function 
renormalization with some details given in the appendix. As well as in the case of electron vertex 
corrections, the nucleon vertex will have an infrared divergence at the pole q — > (for the proton). 
A detailed discussion of the treatment of this type of divergence is given in jsj]. 



III. NUMERICAL RESULTS AND CONCLUSIONS 



The application of the methods described above for {e — N} scattering can be found in calcula- 
tions of the weak charges of the nuclei. Consider the parity- violating Hamiltonian in Eq.j9]). Here, 
for a heavy nucleus we have a coherent effect for V (N) (g) A (e) : 

(u N ^u N ) -> p nuc (r) 5^0. (22) 

The contribution coming from V (e) <g> A (N) is small as it depends on unpaired valence nucleons. 
The latter determines the Hamiltonian for the electron parity-violating interaction with the nucleus 
in the following form: 

G 

H (r) = -^j^Q W eakl<oPnuc {r) ■ (23) 

A relation between the weak charge Q wea k and form factors {Ci p ,Ci n } is straightforward: 

Qleak = 2C lp (Q 2 ^0GeV 2 ] 



Qleak = 2C ln (Q 2 ^0GeV 2 



(24) 



If we take into account only the leading order of the interaction, the weak charge of the proton and 
neutron have the simple definitions: 



O m , = I -As 2 

^weak x ^wi 



n n(0) _ 1 
^weak " - 1 ' 



(25) 



and for the nucleus 



Q W eak = Z-Q p weak + N-Ql eak , (26) 

where Q p wea },, Qweak are weak charges of the proton and neutron including NLO corrections. 

As one can see from Fig- ((TJ) , the input of the NLO corrections is much more significant for the 
proton than for the neutron. This makes experiments involving neutrons more interesting, as they 
allow for the cleaner weak charge extraction. 

Also, form factors C\ v and C\ n normalized by the monopole term jr^m (see Fig.©) do not 
exhibit strong momentum transfer dependence and hence are model independent. This justifies to 
some extent our choice of the monopole form factor approximation. 

Extrapolation of the C\ v and C\ n to zero momentum transfer point (see Fig.©) gives: C\ v = 
0.0360 ± 0.0005 and C Xn = -0.4950 ± 0.0020. The most important of these numbers is the result 
for C\ p . It has been evaluated most recently in [2j, leading to C\ v = 0.0358 db 0.0003. 
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Figure 7: Momentum transfer dependence of the form factors {Ci p , C\ n } and {Ci p , Ci n } ■ A q-P . The 
dotted line represents parity-violating form factor at the tree-level (LO) only, and the solid line shows the 
result with NLO contribution included (LO+NLO). 



The available results from atomic parity-violating experiments for the weak charges of 
Cslf, TZ|? 5 and Bijf ( J3] and references within, and 0) can be used as an experimental test 
of the theoretical predictions: 



C74f(exp) 
T/ 8 2 ? 5 (exp) 



-72.65 ±0.28 ±0.34, 
-114.8 ± 1.2 ±3.4, 



(27) 



Bigf(exp) = -140 ±40. 

Here, the errors are statistical, systematic and coming from an uncertainty of the atomic-physics the- 
ory, respectively. For example, in the case of Cs^ 3 , one should observe 7s (excited) — > 6s (ground) 
parity- violating electric dipole transitions in order to extract the weak charge of Cslf, 3 . This requires 
an accurate knowledge of the atomic wave functions. Some of the analysis of the issues involved, 
with an extensive reference list, is given in [l6| and Using Eq. (l2"4"]) and Eq. ([26l) we compute the 
following results for the corresponding nuclear weak charges: 

Csi 3 (theory) = -73.26 ± 0.26, 

(28) 

T/fj' 5 (theory) = -114.9 ±0.4, 



Bi^f (theory) = -118.8 ±0.4, 

which clearly agree with the experiment. The theoretical uncertainty is estimated quite generously 
and comes mostly from the numerical integration and the extrapolation to the zero momentum 
transfer point. Although the NLO corrections contribute only ~ 1.5% to the results in Eq. PSj) . 
their value is still about four times larger than our theoretical uncertainty. Thus, the more precise 
is the experiment (like Cslf 3 ), the more important it becomes to evaluate most carefully the full set 
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of the applicable NLO corrections. Although our model predictions for the nuclear weak charges 
are in good agreement with the available experimental results, to allow more definitive conclusions 
about the validity of the proposed computational model, the experimental errors for precision 
measurements of the weak charge of the nucleon will have to be reduced. 

Any significant deviation of the weak charge of the proton from the Standard Model prediction 
at low Q 2 would be a signal of new physics. The proton's weak charge Q v wea j, is a well-defined 
experimental observable. At Q 2 — ► the asymmetry can be parametrized as 

[Q 2 Q p weak + Q 4 B(Q 2 )}, 

where B{Q 2 ) is a function of Sachs electromagnetic form factors M related to the Dirac and 
Pauli form factors by the following expression: 

Gl = Fl ~^ F2 ' G m = F ^ + F2 - (M) 

The measurement of Q p weak will be done by the Qweak collaboration fl7l |. and may have extremely 
interesting physical implications. See, for example, [2J]. It will allow one to determine the proton's 
weak charge with ~4% combined statistical and systematic errors, which leads to 9a for the running 
of the weak mixing angle. The Standard Model evolution predicts a shift of /S.sin 2 9w = +0.007 at 
low Q 2 with respect to the Z° pole best fit value of 0.23113 ±0.00015. The weak mixing angle at the 
energy scale close to the Z° pole was measured very precisely, but a precision experimental study of 
the evolution of sin 2 6w to lower energies still has to be carried out. The asymmetry measurements 
proposed for Q wea k experiment will go as low as Q 2 = 0.03 GeV 2 , making it a very competitive 
experiment. 

Using the Dirac and Pauli form factors, we can now compute the extensive set of one-loop 
hadronic electroweak NLO corrections along with the weak charges of the proton and neutron. 
Using a monopole approximation for the form factors, we modified general electroweak couplings by 
inserting the appropriate form factors into the vertices. The monopole structure of the form factors 
allows us to substitute one Feynman diagram with "structured" nucleon with a set of diagrams 
involving only point-like nucleon. This expansion can be visualized by adding fictitious additional 
vector bosons of mass A. For renormalization, we choose the on-shell renormalization scheme. 

In conclusion, it is evident that computational symbolic packages such as FeynArts and FormCalc 
can be efficiently adapted for the theoretical evaluation of NLO effects in electron- nucleon scattering. 
Refs [3| and [lii] give some additional details. A newly proposed measurement of the electron weak 
charge in parity-violating Moller scattering at 12 GeV at JLab could be a very interesting challenge, 
for example. 
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IV. APPENDIX 



A. Some analytical details on couplings 

To introduce couplings between vector bosons and nucleon into the model files of FeynArts, 
instead of vector and axial-vector representation, we employ use of chirality projectors. That can 
be achieved if we compare Eq.JT]), Eq.([2]) with Eq.JU), combined with Eq.© and Eq.JT}, then it is 
possible to write 



(q) = ie[ 



9r (g)+gr v (g) ., , 9jr ( g ) - at- »_ (g) 

n 7|t "T 



7^75 + 



where 



+ 



2m, 



<V*g a /2 (g)] 



IT" (g) = ie 



9t N (g) 7^- + 9 Z R N (g) 7^+ + ^"mth (g) 

gZ' jV (g) + ^' jV (g) 



5l N (g) 7/^- + 5fl JV (g) 1^+ + ir—v^q a F2 (q) 



7-iV 



2rriN 

-c 

I 

i 



^ + 2^^ 



2rriN 



9 f;/(g) = ^-(Fr (iY) (0)-4 S 2 F l( 0)±^(0) 



A 2 



A 2 - g 2 



22/ (g) = *T " (g) = F 1 (0) 



A 2 



A 2 - o 2 ' 



(30) 



(31) 



(32) 



In Eqs. (l30j) and JHTJ, we have adopted the general structure of the coupling from Eq.JT]), with 
coupling strengths g^~ R N (q) and g]^ R (q) given by Eq. ([32l) , Vector part of the coupling a^ a q a also 
can be replaced by the chirality projectors oj± in the following way 



0>ag 



{[-fa, A] ™- + [lu,,A] 



(33) 



Moreover, to adopt FeynArts representation of the coupling through the product of generic and 
class type of the couplings, we introduce the following matrix representation of Eq. ([30]) and Eq. fl3T|) : 

r {N, N , V/j,) = y 7^ro_ , 7 M ro+, [7 M , A] ™- , [in, A] ™+ ) (*NNV (g) , (34) 

with ~C!nnv defines coupling between classes of the vector bosons and nucleons and expressed as a 
2x4 matrix 



G 



nnv = ie 



( gl-(q) GXl N \ 

9 V n N (g) GXK N 
-^Ft N ( q ) 



(35) 



13 



The second column of G nnv represents counterterms of the first order. The Pauli form factor 
pV-N - n Eg. (f35|) has the following structure: 



E 



Z-N 



(?) = h(q) 



E 



rf-N 



(36) 

2 (?) = Fi (q) ■ 

The coupling defined in Eq. (|35|) has a counterterm part at the one-loop level represented by the 
matrix, which has following structure: 

' r Re[Sfu 
r Re[5f R ] 

(q) Re[5f L ] 



( gXe n \ 


( 


n V-N 
^lR 




n V-N 
2L 




\ r V-N j 
\ Lr 2R I 


\ 



9 V L ~ N 
9 V R~ N 

1 jpV-N 
' im N r 2 

V-N 



1 TP 

4rojv 2 



(37) 



(q)Re[5f R )J 



where the hadronic field renormalization constants 5 /l,r are computed using the expansion on 
Fig ®- More details are given in Ref.| 





Figure 8: Hadronic wave function renormalization expansion. 
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